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Abstract. Diffusive shock acceleration (DSA) at relativistic shocks is expected to be an important 
acceleration mechanism in a variety of astrophysical objects including extragalactic jets in active 
galactic nuclei and gamma ray bursts. These sources remain strong and interesting candidate sites for 
the generation of ultra-high energy cosmic rays. In this paper, key predictions of DSA at relativistic 
shocks that are salient to the issue of cosmic ray ion and electron production are outlined. Results 
from a Monte Carlo simulation of such diffusive acceleration in test-particle, relativistic, oblique, 
MHD shocks are presented. Simulation output is described for both large angle and small angle 
scattering scenarios, and a variety of shock obliquities including superluminal regimes when the 
de Hoffman-Teller frame does not exist. The distribution function power-law indices compare 
favorably with results from other techniques. They are found to depend sensitively on the mean 
magnetic field orientation in the shock, and the nature of MHD turbulence that propagates along 
fields in shock environs. An interesting regime of flat spectrum generation is addressed, providing 
evidence for its origin being due to shock drift acceleration. The impact of these theoretical results 
on gamma-ray burst and blazar science is outlined. Specifically, Fermi gamma-ray observations 
of these cosmic sources are already providing significant constraints on important environmental 
quantities for relativistic shocks, namely the frequency of scattering and the level of field turbulence. 

Keywords: Diffusive shock acceleration; hydromagnetic turbulence; gamma-ray bursts, active 
galactic nuclei, non-thermal emission 

PACS: 98.70.Rz; 95.85.Pw; 98.70.Sa; 52.35.Ra; 52.25.Xz; 52.27 .Ny; 52.35.Tc; 52.65.Pp 

INTRODUCTION 

There is bountiful evidence for efficient particle acceleration at collisionless shocks in 
the universe. The heliosphere, with its planetary bow shocks and traveling interplanetary 
shocks, has provided interesting and useful test cases for shock acceleration theories. In 
the remote regions of the universe, supersonic jets from active galactic nuclei (AGNs) 
and gamma-ray bursts (GRBs) have offered fascinating windows into an energetic part 
of the cosmos, where non-thermal radio waves, X-rays and gamma-rays abound. Fully 
understanding these sources mandates knowledge of how the particles that generate their 
light emission are energized. The foremost paradigms invoke acceleration at relativistic 
shocks. Radio imaging reveals very structured and time-variable jets in AGNs, rapid X- 
ray and gamma-ray variability in both blazar AGNs and GRBs suggest compact emission 
regions with relativistic bulk motions. Accordingly, comprehending the relationship 
between shock acceleration predictions and observations of these sources offers the key 
to elucidating the understanding of their environs; this constitutes the focus of this paper. 



To effect this, an investigation of the features of diffusive shock acceleration is 
presented, using results from a test particle Monte Carlo simulation. Then the paper 
addresses probes of the acceleration theory parameter space imposed by extant GRB and 
blazar observations in high energy gamma-rays. The Monte Carlo approach JUS 0,111] 
is one of several major techniques devised to model particle acceleration at relativistic 
shocks; others include semi-analytic solutions of the diffusion-convection equation [0,0, 
8], and particle-in-cell (PIC) full plasma simulations [@, 1(1 11. 12]. Each has its merits 
and limitations. Tractability of the analytic approaches generally restricts their solutions 
to power-law regimes for the phase space distributions /(p) . PIC codes are rich in their 
information on shock-layer electrodynamics and turbulence. However, to interface with 
astrophysical spectral data, a broad dynamic range in momenta is desirable, and this is 
the natural niche of Monte Carlo simulation techniques. A core property of acceleration 
at the relativistic shocks is that the distribution functions /(p) are inherently anisotropic. 
This renders the power-law indices and other distribution characteristics sensitive to 
directional influences, such as the magnetic field orientation with respect to the shock 
normal, and the nature of MHD turbulence that often propagates along the field lines. 
These connections between observables and physical parameters of the shock environs 
are studied in some depth in this exposition. 



DIFFUSIVE ACCELERATION AT RELATIVISTIC SHOCKS 

The exploration of the properties of diffusive shock acceleration that is presented here 
employs the well-known kinematic Monte Carlo technique of Ellison and Jones that 
has been employed extensively in supernova remnant and heliospheric contexts, and 
is described in detail in numerous papers [Q], 0, [3 EH EBL Q2]. Particles are injected 
upstream and allowed to convect into the shock, meanwhile diffusing in space so as to 
effect multiple shock crossings, and thereby gain energy through the shock drift and 
Fermi processes. In general, the upstream fluid frame magnetic field is inclined at an 
angle Bfl to the shock normal. The particles gyrate in laminar electromagnetic fields, 
with their trajectories being obtained by solving the Lorentz force equation in the normal 
incidence shock rest frame (NIF), where the upstream flow is incident along the shock 
normal (x-direction: see the left hand panel of Fig.[[]for the NIF shock geometry). In this 
frame, when Bfl > , there is a u x B drift electric field in addition to the magnetic 
field. The effects of Alfven wave and other hydromagnetic magnetic turbulence on 
particle propagation are modeled by phenomenologically scattering these ions elastically 
in the rest frame of the local fluid flow. This is generally applicable for high Alfenic 
Mach number shocks. The simulation outputs particle fluxes and momentum and angular 
distributions, usually in the NIF, at any location upstream or downstream of the shock. 

The simulation can routinely model diffusion incurred in particle interactions with 
MHD turbulence using either large-angle or small-angle scattering. At every scattering, 
the direction of the particle's momentum vector p is deflected in the local fluid frame 
to a new value p'. The scattering angle 80 = cos _1 (p.p'/|p| |p'|) is uniformly sampled 
within a solid angle 8£l up to a maximum deflection angle scatt (see Figs. 2 and 3 
of [0] for the scattering geometry). The time, 8tf, between scatterings in this frame is 
coupled yj] to the mean free path, A, and scatt , via 8tf ps A0 S c a tt/( ov ) f° r particles 




FIGURE 1. The geometry in the normal incidence (NIF; left panel) and de-Hoffmann Teller (HT; 
right panel) shock rest frames. Upstream flow speeds in the two reference frames are related by 
miht = )3ihtc = Mi/cos©Bfi- Upstream and downstream quantities are denoted by subscripts 1 and 
2 , respectively. In general, the NIF field angle ©b s i differs from the fluid frame/HT frame value ©Bfi ■ 
Also, the shock plane in the HT frame is rotated from that in the NIF due to relativistic aberration effects. 



of speed v . The resulting effect is that the gyrocenter of a particle with gyroradius r g 
is shifted randomly by a distance of the order of scat t^g in the plane orthogonal to the 
local field. Accordingly, cross-field diffusion emerges naturally from the simulation. For 
large angle scattering (LAS, defined more precisely below), the scattering solid angle 
is sca tt < ^ steradians, and the transport is governed by kinetic theory JTsLCLtI], where 
the ratio of the spatial diffusion coefficients parallel (Mi = Av/3) and perpendicular 

( K± ) to the mean magnetic field is given by K±/Kn = 1/(1 + r\ 2 ) . Here, the parameter 
r\ = X/r g is the ratio of a particle's mean free path X to its gyroradius r g . As is 
often implemented for simplicity, in this work X is assumed to be proportional to 
the particle momentum p, so that r\ is independent of p. This can be adjusted to 
accommodate other scattering laws as desired. Clearly, r\ controls the amount of cross- 
field diffusion, and is a measure of the level of turbulence present in the system, i.e. is 
an indicator of (8B/B) . The Bohm limit of quasi-isotropic diffusion is realized when 
r\ ~ 1 and (8B/B) ~ 1 . This phenomenological description of diffusion in Monte Carlo 
techniques is most appropriate at high energies, and omits the details of microphysics 
present in plasma simulations such as PIC codes. In the injection domain at slightly 
suprathermal energies, the influences of complex turbulent and coherent electrodynamic 
effects become important, and will substantially modify the picture from pure diffusion. 

It will become clear below that in oblique relativistic shocks, the diffusive transport of 
particles across the field, and hence across the shock, becomes critical to their retention 
in the acceleration process. Accordingly, for such systems, the interplay between the 
field angle Bfl and the value of X /r g controls the spectral index of the particle distri- 
bution |0,@], a feature that is central to the interpretation of astrophysical source spectra. 
The phase space for acceleration theory results is now explored, before addressing their 
relevance to two classes of cosmic sources: gamma-ray bursts and blazars. 



ACCELERATION SIGNATURES AT RELATIVISTIC SHOCKS 



Diffusive particle acceleration at relativistic shocks naturally divides into two regimes: 
mildly-relativistic (u\ x < c) and ultra-relativistic (u\ x ~ c) systems. Hereafter, sub- 
scripts 1 denote upstream quantities. Neither regime is accessible to diagnostics via in 
situ particle measurements since the principal astrophysical sources possessing relativis- 
tic shocks exhibiting signatures of acceleration, namely pulsar wind termination shocks, 
jets in active galactic nuclei (blazars) and gamma-ray bursts, are so distant. An impor- 
tant property of diffusive acceleration at relativistic shocks that distinguishes them from 
their non-relativistic counterparts is their intrinsic anisotropy. This is driven by the pow- 
erful convective influence that enables efficient loss of particles away and downstream 
from the shock. The immediate consequence of this loss is a general difficulty in gener- 
ating flat distributions of shock-accelerated particles, particularly for so-called superlu- 
minal (/3i x / cos0 Bfl > 1 ) relativistic discontinuities. These losses are muted somewhat 
in mildly-relativistic shocks, which are naturally expected in GRBs and perhaps also in 
blazars. Moreover, the resulting distributions dN /dp p 2 f(p) are more or less com- 
mensurate with those required to explain the radiation observations from these sources. 
Accordingly, this paper will study T\ = (1 — u 2 x /c 2 )~ 1 / 2 ~ 1 — 5 shocks, focusing on 
the power-law tails dN /dp p a realized at high momenta, where all memory of the 
injection momentum has been destroyed by diffusion. The array of possible distribution 
indices o is highlighted, spawned by the sensitivity of both the energization in, and es- 
cape from, the shock layer, to (i) the size of the momentum deflection angle <9 scatt , (ii) 
the frequency or relative mean free path X /r g of scatterings, and (iii) the upstream field 
obliquity Bfl , a quantity connected to the global MHD structure of the shock. 

The first effect, the dependence of o on scatt , has been documented in several recent 
papers [0, 0, @], and will just be summarized here. When the diffusion in the shock 
layer samples large field fluctuations 8B/B ~ 1 (seen, for example, in PIC simulations 
of relativistic shocks driven by the Weibel instability |@, Oil CLll it corresponds 
to large momentum deflections, delineating the regime of large angle scattering (LAS) 
with A/Ti < scatt < n), where Ti is the upstream flow's incoming Lorentz factor. This 
regime was first explored for Ti < 5 by Ellison et al. yj]. Such large deflections produce 
huge gains in particle energy, of the order of T 2 , in successive shock crossings. These 
gains are kinematic in origin, and are akin to those in inverse Compton scattering. The 
result is an acceleration distribution dN /dp that is highly structured and much flatter on 
average yj] than p~ 2 for strong, parallel (0 B n = 0°), shocks, i.e. those with high sonic 
and Alfenic Mach numbers. The bumpy structure is kinematic in origin, corresponding 
to sequential shock transits lfl8ll . and becomes more pronounced I^MEl [l9ll for large 
Ti . For ultra-relativistic shocks, information of the injection momentum scale becomes 
insignificant when p ^> mc, and the bumps asymptotically relax to form a power-law 
distribution dN /dp p~° , with an index in the range of o ~ 1.6 [0]. From the plasma 
physics perspective, magnetic turbulence in relativistic shocks could easily be sufficient 
to effect scatterings on intermediate to large angular scales scatt > 1 /Yi , a proposition 
that becomes more enticing for ultra-relativistic shocks. 

The principal focus in this paper is on particle distributions for scatt < 1 /Y\ regimes, 
which are much smoother in appearance, and often necessarily steeper, at least for su- 



perluminal regimes. The property of distribution smoothness meshes more easily with 
radiation spectral observations of extragalactic astrophysical sources, thereby motivating 
exploration of this portion of phase space. The kinematic energy gains in shock cross- 
ings are lowered considerably [|18|l when scatt drops below this "Lorentz cone" angle 
1/Ti. Accordingly, the character of dN/dp and the particle anisotropy at the shock 
dichotomize, partitioned by the scatt ~ 1 /T_\ boundary. Intermediate scattering angles 
S catt ~ 1 /Ti generate smooth distributions lULQil], much like those for small angle scat- 
tering (SAS, often called pitch angle diffusion, PAD). The SAS regime has spawned the 
often cited asymptotic, ultrarelativistic index of o = 2.23 for dN/dp p a [8], first 
noticed in Monte Carlo simulations Il20ll . This special result, applicable in both shock 
rest and fluid frames for momenta p 3> Time, is realized only for parallel shocks with 
®Bfi = 0° in the limit of scatt <C 1 /T\ , where the particle momentum is stochastically 
deflected on arbitrarily small angular (and therefore temporal) scales. In such cases, 
particles diffuse in the region upstream of the shock only until their velocity's angle 
to the shock normal exceeds around 1/Ti , after which they are rapidly swept down- 
stream of the shock. The lower kinematic energy gains in shock transits dominate higher 
shock-layer retention rates, and guarantee a steeper distribution under SAS conditions 
for Bfl = 0° shocks; the monotonic steepening and loss of structure as scatt declines is 
exhibited in Q4f| for ultra-relativistic shocks, and Ill9ll for the mildly-relativistic domain. 

Now the focus turns to the influence the effective frequency X/r g of scatterings, 
and the upstream field obliquity Bfl have on the accelerated population. Represen- 
tative particle (electrons or ion) differential distributions dN/dp that result from the 
simulation of diffusive acceleration at mildly-relativistic shocks of speed /3i x = 0.5 
are depicted in the left panel of Figure |2] (see |0, 3] for T\ ^> 1 simulation results). 
Here, the subscript x denotes components along the shock normal. These distributions 
were generated for scatt < 10°, i.e. in the SAS regime, for low magnetic fields corre- 
sponding to Alfvenic Mach numbers M A ^> 1 , and in the NIF frame. Results are dis- 
played for two different upstream fluid frame field obliquities, namely B fi = 48.2° 
and Bfl = 59.1°, with corresponding de Hoffman-Teller (HT 112 ill ) frame dimension- 
less speeds of j8 1H t = Pix/ cos0 Bfl = 0.75 and 0.975, respectively. Subluminal shocks 
are those where the HT flow speed /3 1H t corresponds to a physical speed, less than 
unity, i.e. the upstream field obliquity satisfies cos0 Bf i < j5\ x . When j8i H t > 1 » the de 
Hoffman-Teller frame does not exist, and the shock is said to be superluminal. See the 
right hand panel of Fig. Q] for the HT frame shock geometry. The distributions clearly 
exhibit an array of indices o , including very flat power-laws, that are not monotonic 
functions of either the field obliquity Bfl nor the key diffusion parameter r\ = X/r g . 
These properties are illustrated in the right panel of Fig. |2l where it is also evident that 
the distributions are generally steeper in superluminal shocks [2] with A/r g > 10. The 
left panel of Fig. [2] also emphasizes that the normalization of the power-laws relative to 
the low momentum thermal populations is a strongly-declining function of X/r g . This 
is a consequence of a more prolific convection of suprathermal particles downstream 
of the shock that suppresses diffusive injection from thermal energies into the acceler- 
ation process. Such losses are even more pronounced when X/r g > 10 4 , to the point 
that acceleration is not statistically discernible for j8 1H t > 0.98 runs with 10 4 simulated 
particles. This feature is salient for the discussion on astrophysical sources below. 




FIGURE 2. Left panel: Particle distribution functions dN /dp (scaled by p 2 ) from mildly-relativistic 
sub-luminal shocks (T\ x fi\ x = 0.577, i.e. /3u = U\ x /c = 0.5) of upstream-to-downstream velocity com- 
pression ratio r = u\ x J ii2 x ~ 4 . Simulation results are depicted for two de Hoffman-Teller frame upstream 
flow speeds /3iht — Pix/ cos ©Bfi ■ These are in distinct groups of four: /3iht = 0.75 (upstream magnetic 
field obliquity ©Bfi = 48.2° , multiplied by 10 3 ) for the upper four histograms, and /3iht = 0.975 (i.e., 
©Bfi = 59.1°) for the lower four histograms. Scattering off hydromagnetic turbulence was modeled by 
randomly deflecting particle momenta by an angle within a cone, of half-angle sca tt , whose axis coin- 
cides with the particle momentum prior to scattering; four different ratios of the diffusive mean free path 
X to the gyroradius r g were adopted for each ©Bfi ■ AH results were for small angle scattering (SAS), 
when sca tt < 1 /Ti and the distributions become independent of the choice of sca tt ■ 
Right panel: Power-law indices a for simulation runs in the limit of small angle scattering (pitch angle 
diffusion), for mildly-relativistic shocks of upstream flow speed $\ x = U\ x jc = 0.5 , and an MHD veloc- 
ity compression ratio r = 4 . The indices are displayed as functions of the effective de Hoffman-Teller 
frame upstream flow speed /3iht = Pu/ cos ©Bfi > with select values of the fluid frame field obliquity 
©Bfi marked at the top of the panel. The displayed simulation index results were obtained for different 
diffusive mean free paths X parallel to the mean field direction, namely X jr g = 1 (squares), X jr g = 10 
(triangles), Xjr g — 10 2 (pentagons), and X/r g = 10 3 (triangles), as labelled. The lightweight curve at 
the bottom labelled KH89 defines the semi-analytic result from Kirk & Heavens' |7] solution to the 
diffusion-convection equation, corresponding to X jr g — > °° . The short heavyweight lines labelled GRB 
930131 (EGRET detection) and GRB 080916c (Fermi detection) indicate the approximate spectral index 
a that is appropriate for these gamma-ray bursts, if a cooled synchrotron emission scenario is operable. 



The existence of very flat distributions in the subluminal domain for very large 
A /r g is a remarkable feature of Figure [21 This phenomenon was identified by Kirk & 
Heavens [7] in their eigenfunction solution technique [6] for the diffusion-convection 
equation, which was restricted to subluminal, oblique shocks. Results for their /3i x = 0.5 
analysis are presented as the solid curve labelled KH89 in the right panel of the Figure. 
Clearly, the Monte Carlo indices closely approach those of the semi-analytic method of 
[0] for X/r g = 10 3 , an agreement that is improved slightly when A/r g is increased 
to 10 4 . Yet there are differences between the two approaches, and these yield the 
expected slight discrepancies in spectral index determination. The method of Q7Q, being 
tantamount to a guiding center technique, employs conservation of the magnetic moment 



p 2 (l — /i 2 )/ (25) for particle-shock interactions when determining transmission and 
reflection probabilities of charges, naturally differing from gyro-orbit determinations of 
these probabilities. This nuance probably seeds the different values of o derived in the 
two approaches, which, as expected, are small when shocks are effectively parallel, and 
highly oblique. More details of the comparison of these approaches are discussed in Rial , 
which provides an extended exposition on the implementation of the Monte Carlo code, 
its validation and many of its key acceleration results. It can also be noted that comparing 
with the indices derived by yVD was the major motivation behind the artificial choice 
of the compression ratio r = 4 , which is somewhat larger than the Rankine-Hugoniot 
MHD value for /3i x = 0.5, M s = 4.04 conditions. Also, a low sonic Mach number M s 
was chosen so as to maximize the efficiency of injection from thermal energies. 



The Action of Shock Drift Acceleration 

The origin of these flat indices is elucidated by the parameter survey here, and more 
specifically by the inclusion of cross field diffusion in the Monte Carlo simulations. 
Diffusive transport of particles perpendicular to the mean field was omitted in the 
analysis of [7], a restriction that is the primary reason for o dropping below two. This is 
clearly evident from the right panel of Figure [2l which indicates that a < 1.5 scenarios 
are realized only for A/r g > 10 2 , i.e. laminar fields and almost pure gyrational motion. 
The primary origin of the acceleration is then connected to coherence in the shock layer. 
To provide further insight, individual particle trajectories were tracked in the Monte 
Carlo runs, and those exhibiting profound energy gains isolated. One such example is 
provided in Fig. [31 together with its corresponding momentum trace (right panel); the 
particles were injected with superthermal momenta p = 5mc to circumvent any injection 
problems. Diffusion has minimal impact on the gyrational motion in the left panel. The 
trajectories are 2D projections, and the pitch angle is evolving almost adiabatically in the 
particle-shock interaction to preserve gyrational coherence. The momentum histories in 
the right hand panel more readily highlight the property that the acceleration is directly 
coupled to periods when the particle's gyration straddles the shock. Moreover, inertial 
motion in the y or uxB direction accompanies these epochs of energization. These two 
characteristics are the hallmarks of shock drift acceleration (SDA). 

The acceleration of particles in the drift electric fields associated with oblique shock 



discontinuities has been extensively studied in non-relativistic contexts [|22L I23L 1241 . 125 
HfJ. The origin of the effect is the net work done on a charge by the Lorentz force in a 
zone of non-uniform magnetic field. The principal equation governing this is 



P # = 4P-( E+ c xB j = 4P' E > (!) 

where E is the uxB drift field that exists in any oblique shock rest frame other than the 
HT frame. In the uniform B fields either upstream or downstream, the energy gains and 
losses acquired during a gyroperiod exactly cancel, so that no net work is done, dW = . 
In contrast, when a charge's gyromotion straddles the shock discontinuity, the sharp field 
gradient induces an asymmetry in the time spent by the charge either side of the shock, 




FIGURE 3. Left panel: A portion of a shock-layer trajectory for a selected particle accelerated from low 
to high energies, illustrating two key features of shock drift acceleration: coherent trapping in the shock, 
and occasional upstream excursions. The projection is onto the x-y plane, where the u x B drifts lie in 
the y-direction. The shock speeds are as in Fig. [2 the large value of X/r s — 10 4 is adopted to generate 
the requisite conditions for coherent interaction with the shock. In this projection, energy gains incurred 
while the particle's gyrations straddle the shock only become obvious when its pitch angles are altered in 
upstream transits. Particles that undergo these trajectories are rare, but populate the power-law tail. 
Right panel: The coupling between the fluid frame (p F ) and NIF shock rest frame (ps) momenta of 
the particle selected in the left panel, and the drift coordinate y. The shock frame momentum evinces 
oscillatory behavior during upstream excursions, and both p F and p s display the linear trend with drift 
y that is characteristic of shock drift acceleration. The fluid frame momentum exhibits a "rectangular 
hysteresis" while the particle gyrates in the shock layer. 



so that energy gains and losses do not negate each other. The compressive nature of the 
field discontinuity biases the net work done to positive increments in shock encounters 
between upstream excursions, and it is simply shown 11231 . 12511 that dW = qE y dy, i.e. this 
energy gain scales linearly with displacement along the drift coordinate y . This is the 
punchline of the right panel of Fig. |3l where y effectively represents a time coordinate 
during shock drift episodes, with dW /dt = qE y / c . It is also noteworthy that interspersed 
between these acceleration periods are upstream excursions where infrequent scattering 
slowly tries to isotropize the pitch angles, a feature clearly identified for the SDA phe- 
nomenon at non-relativistic shocks Il26ll . Particles that participate in the SDA initially 
have HT frame pitch angle cosines /i H t considerably less than unity, so as to satisfy the 
reflection criterion: therefore they are rare in near-luminal shocks, for which the inci- 
dent upstream angular distribution is highly-beamed around jU H t ~ 1 • During upstream 
excursions that follow reflection, the particles tend to get swept back to the shock as 
soon as they are deflected and acquire momenta outside their Lorentz cone, i.e. before 
reaching isotropy in the upstream fluid frame. Therefore, in these X/r g ^> 1 cases, jU HT 
continues to satisfy the reflection criterion during these brief upstream epochs. The rep- 
etition of SDA and upstream excursions, i.e. trapping in and near the shock, is therefore 
virtually guaranteed once the initial reflection is realized for a select particle. 



This investigation provides an identification of the significance of shock drift accel- 
eration in controlling o for relativistic shocks. The reason it couples to unusually flat 
distributions (o < 1.5) revolves around extremely efficient trapping in the shock layer, 
which permits repeated episodes of SDA in select particles with appropriately-tuned 
gyrophases at the onset of shock-orbit interactions. The combination of the rapid ener- 
gization rate, upstream hiatuses in SDA, and a slow leakage rate downstream leads nat- 
urally Il27n to an approximate dN /dp p~ l distribution. But only when k/r g > 10 3 . 
The introduction of turbulence easily disrupts the coherence and precipitates efficient 
convection downstream [27], quenching the effectiveness of SDA. The result then is 
a dominance of shock drift acceleration contributions by first order Fermi (turbulent) 
ones when A/r g < 30. More importantly, when k/r g > 10, an inexorable sweeping of 
charges downstream in superluminal shocks overpowers both acceleration contributions 
and steepens the spectrum dramatically, as is evident in the right panel of Fig. [2l 



ASTROPHYSICAL SOURCE CONTEXT 

The shock acceleration theory results presented in the previous Section can now be in- 
terpreted in the light of observations of astrophysical sources. The first class of germane 
sources consists of gamma-ray bursts (GRBs), whose prompt emission is generally ob- 
served in the 10 keV - 10 GeV range, with a characteristic spectral break around 200 
keV to 1 MeV. Above this break, the spectrum is generally (but not always) an extended 
power-law driy/d£y £y a ' ! . The focus here is on the relationship between the high en- 
ergy spectral index and the underlying particle acceleration conditions. This forges 
a direct connection to data from CGRO's EGRET telescope, and now to the growing 
database of Fermi LAT burst detections. The EGRET Of/, index distribution |28|] is con- 
stituted by a handful of sources with indices scattered in the range 2 < < 3.7 , with 



brighter bursts' indices concentrated in the range Of/, < 2.8, as tabulated in [29]. The 
recent Fermi detection S of GRB 080916c in both the GBM and LAT instruments 
offered an index of «/, ~ 2.2 at energies above ~ 2 MeV in its most luminous epoch, 
and a steeper spectrum ( a% > 2.5 ) at other times. It is then evident that observationally, 
shock acceleration models must accommodate a radiation spectral index in the range 
2 < ah < 4 in order to be viable. Moreover, they must reasonably account for the spec- 
tral variability identified in GRB 080916c, i.e. fluctuating values in a given source. 

If one presumes that these photon spectra result from synchrotron emission that 
rapidly cools the radiating electrons, a popular paradigm yU, [32l Q3[] for the production 
of the prompt emission, then ah = (c + 2) /2 for electron acceleration populations 
dN I dp p ° . The o corresponding to the values of «/, = 2.0 for the flat-spectrum 
EGRET burst GRB 930131 and a h = 2.2 for GRB 080916c for time intervals (c,e) 
| f30h are marked on the right panel of Fig. [2] It then becomes clear that the emission 
in GRB 930131 is consistent with acceleration at subluminal shocks and relatively near 
the Bohm diffusion limit. In contrast, the slightly steeper GRB 080916c spectrum is 
better explained by mildly superluminal shocks in the 60° -70° obliquity range, but 
only if the scattering is strong, i.e. A/r g < 3 — 10. If, instead, synchrotron cooling is 
inefficient, the photon differential spectral index given by ah = (o + l)/2. Then the 
injected distribution must have an index a higher by unity than that for cooling models, 



in order to match the burst observations. This is a profound difference in that it moves the 
viable shock parameter space into the superluminal range, i.e. at higher field obliquities, 
and Bohm-limited diffusion is observationally excluded. No bursts have so far evinced 
extended power-law spectra flatter than GC/, « 1.95, absolving the need for acceleration 
in shocks with extremely low turbulence, i.e. A/r g > 10 2 regimes. This is fortunate, 
since, from the left panel of Fig. [2l such shocks are inherently inefficient accelerators. 
Moreover, the generation of field turbulence is a natural part of dissipation in shocks, 
so that almost laminar fields are not expected, nor observed in in situ magnetometer 
measurements at heliospheric shocks (e.g. rf34l l3~5ll and references therein). 

The second relevant astrophysical context concerns blazars, the subset of active galac- 
tic nuclei possessing relativistic jets of material emanating from the supermassive black 
holes at their centers; these jets are oriented virtually towards the observer. These have 
been the preserve of gamma-ray experiments ever since their discovery by EGRET 
ll36ll . and subsequent observation by ground-based Cherenkov telescopes at TeV en- 
ergi es |37|1 . The TeV-band signals typically exhibit steep photon spectra (e.g. Mkn 421: 
B8U 39II') that include the absorption due to pair producing interactions 77^ e + e~ with 
infra-red and optical light generated by the intergalactic medium along the line of sight 
to the observer. The correction for this attenuation (so called de-absorption), leads to 
the inference of extremely flat particle distributions in energetic gamma-ray blazars (see 
for example |4j]), with indices as low as o ~ 1.0 — 1.5 in high redshift sources. To ac- 
commodate these observational constraints would require subluminal shocks with very 
modest or low turbulence levels. Yet, such inferences are based upon measurements in 
a limited waveband subject to profound absorption. More improved diagnostics are now 
enabled by Fermi LAT detections of blazars, which extend the observational window 
over a much larger energy range, nominally from 100 MeV to over 1 TeV, and most 
crucially, below the attenuation window. Accordingly, Fermi observations can more di- 
rectly probe the underlying radiating particle population. A prime example of this is 
the multi- wavelength campaign on the PKS 2155-304 blazar [|40(] in a non-flare state, 
whose combined Fermi-HESS spectrum from 300 MeV to 3 TeV indicates an unatten- 
uated photon spectral index of ot/, ~ 1.6, steepening to «/, ~ 2.0 above 1 GeV For a 
standard inverse Compton scattering interpretation with insignificant radiational cool- 
ing, this translates to an electron power-law index in the range o ~ 2.2 — > 3.0, so that 
acceleration at mildly superluminal oblique shocks should provide the best description. 



CONCLUSIONS 

This paper has investigated some of the key characteristics of particle acceleration at 
relativistic shocks, including the identification of the role of shock drift acceleration 
in generating flat distributions in mildly-relativistic, subluminal shocks. It has also ex- 
plored the connection between the acceleration process and high energy emission in 
two classes of astrophysical sources, namely gamma-ray bursts and jets in blazars. The 
simulation results presented clearly highlight the non-universality of the index of ener- 
getic, non-thermal electrons and ions, spawned by the variety of shock obliquities and 
the character of hydromagnetic turbulence in their environs. This non-universality poses 
no problem for modeling GRB or blazar high-energy power-law indices, though obser- 



vations generally constrain the parameter space to subluminal or highly-turbulent and 
modestly superluminal shocks not far from the Bohm diffusion limit. Diffusive accel- 
eration at ultra-relativistic shocks requires 8B/B ~ 1 in order to generate sufficiently 
low o to mesh with observed source photon spectra. It is unclear whether conditions in 
bursts and blazars can support such levels of turbulence at shocks embedded in their rel- 
ativistic outflows, though it should be noted that large field fluctuations naturally emerge 
in PIC simulations of Weibel instability-driven relativistic shocks. It is anticipated that 
the rich prospects for Fermi gamma-ray observations of blazars and GRBs in the next 
few years will enhance our understanding of particle acceleration in their environs. 
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